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Abstract 

We demonstrate that nonlinear decay of obliquely propagating Langmuir waves into Langmuir 
and Alfven waves (L — > L' + A) is possible in a one-dimensional, highly relativistic, streaming, pair 
plasma. Such a plasma may be in the magnetospheres of pulsars. It is shown that the characteristic 
frequency of generated Alfven waves is much less than the frequency of Langmuir waves and may 
be consistent with the observational data on the radio emission of pulsars. 
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I. INTRODUCTION 



It is generally believed that a coherent (brightness temperatures up to ~ 10^" K) radio 
emission of pulsars is generated by plasma instabilities in the one-dimensional, relativistic 
plasma consisting of electron-positron pairs (e.g., l|). Such pairs may be produced by 
7-rays via electromagnetic cascades in the vicinity of the magnetic poles of pulsars that 
are identified with strongly magnetized (~ 10^ — 10^"^ G), rotation-powered neutron stars. 
Created particles very rapidly lose the momentum component transverse to the magnetic 
field and flow away along the open magnetic field lines with Lorenz-factors ranging from 
7min ~ 10 to 7max ~ 10^. A two-strcam instability has been proposed as a mechanism of 
radio emission of pulsars soon after their discovery 2|. Later, using the available models 
of pulsar magnetospheres it was argued that if the plasma outflow is stationary, the two- 
stream instability does not have enough time to be developed before the plasma escapes the 
pulsar magnetosphere (for a review, see j^). It was argued 0, ^ that the process of pair 
creation near the pulsar surface is strongly nonstationary, and the pair plasma gathers into 
separate clouds spaced by / ~ 10® cm along the direction of its outflow. In this case, at the 
distance ^ '^hmin ~ 10*^ cm from the neutron star fast (7 ^ 7max) particles of one plasma 
cloud overtake slow (7 ~ 7min) particles of the preceding cloud, and mutual overlapping 
of the clouds begins. In the overlapping region, there are, in fact, two steams of slow and 
fast particles, and the condition for the development of two-stream instability is created Q]. 
For typical parameters of pair plasma in the pulsar magnetospheres the growth rate of the 
instability is quite sufficient for its development (see below). 

The frequency of Langmuir (L) waves generated at the distance from the neutron star is 
a few ten times higher than the typical frequency of pulsar radio emission (see [3] and below). 
Therefore, the wave frequency has to decrease significantly because of some processes to be 
compatible with the observational data. In this paper, we discuss the process of conversion 
of L waves into low-frequency Alfven [A) waves that can solve the high-frequency problem. 



2 



II. BASIC EQUATIONS 



We consider pair plasma in the external uniform magnetic field Bq directed along the 
z-axis. The behavior of this plasma may be described by the Vlasov and Maxwell equations, 
d d / I \ d 

g-^fa + v,—U + ea (E, + -[v X B],j — /„ = 0, (1) 

where fa is the distribution function for the particles of type a which is normalized so that 
/ fadpz = riQ = J2a^a IS the particle density, B = Bq + B, B is the magnetic field 
perturbation, and j = J2a J ^fadpz is the current density. 

Our study of nonlinear phenomena will be performed in the frame of the weak turbulence 
theory. The particle distribution functions and the current density may be written in the 
perturbation manner, 

f. = jT + ff^ + fL'^ + .-, 3=f^+3''^+i^ + ..., (3) 

here the unperturbed distribution function f^^ is assumed to be stationary and homogenous. 
Besides, we use the infinite magnetic field approximation (for finite magnetic field effects. 



see jsj). In this approximation, the distribution function is one-dimensional in the velocity 
space (v II Bq) and depends only on Vz- Technically, this means that only the z component 
of j*^*-* is nonzero, 

jf = Y.^JvJ^'^dpz. (4) 
In this case from Eqs. ([l]) and ([3]) we have 

|/r + e„(B.A/<..)^0. (5) 

|/«+.4/<-' + e„(E.|-/r^.)^0. (6) 

where the angular brackets denote an average over one pulsation. 

Performing Fourier transformation with respect to temporal and spatial variables, Eq. 
dHD for i = 1,2 yields 



00 - k^V:. 



— e 
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fi'\k) = I dk'EzCk - k')E^{k r\- ;7- , (8) 

UO — KzVz OPz J UJ — kzVz 



where k = (k, cj). The first (fj^^) and second (fj^^) order perturbations describe the linear 
dynamics of plasma and three wave resonant processes, respectively. 

Without lost of generality we assume that the wave vector k lies in the x, 2;-plane. From 
Eqs. ([2]), (jlj) and ([7]) it can be shown that there are the following three fundamental modes 
of pair plasma (e.g., 9(]). The first one is the extraordinary wave with the electric field 
perpendicular to the k. Bo-plane. The other two modes are L waves generated in the devel- 
opment of the two-steam instability and A waves into which L waves may be converted (see 
below). For both these modes their electric fields lie in the k. Bo-plane. 

In the infinite magnetic field approximation when only one component {Ezz) of the linear 
permittivity tensor is nonzero 9|, 



E 



7^(0; - k,Vz) 



(9) 



for rather low wave numbers, 



the dispersion relations can be written in the form 



uji = 2u'^{r')+^kic 
for quasi-longitudinal L waves, and 

2 ; 2 2 1 



2^2 



+ ktc 



2„2 



ktc 



ky 



(10) 



(11) 



(12) 



2a;2(7(l + t;,/c)2)J 

for quasi-transversal A waves, where Fa = f^V^a^ so that / F^dpz = 1, cUp = (47re^rao/me)^/^ 
is the plasma frequency, and (...) = / ...Fadpz- 

If the conditions (fTOj) are fulfilled for A waves the ratio of their electric field components 

is 

(13) 



2ujl{^{l + vz/cf)- 



III. GENERATION OF L WAVES AND THEIR EVOLUTION 



We consider the model where L waves are generated by the two-stream instability that 
develops at the distance ~ 10^ cm from the neutron star surface because of overlap of 
the pair plasma clouds ejected from the pulsar |l^. Below, all calculations are done in the 



plasma frame where the mean vector-velocity of the cloud particles is zero. We assume that 
in the regions where the plasma clouds overlap the densities of the slow and fast particles 
are the same (for development of two-stream instability for arbitrary distribution function 



of plasma particles, see lOj). In this case the Lorentz-factor of the plasma frame in the 
pulsar frame is 7^ ^ (7max7min)^^^ ^ 10^ 

The two-stream instability starts developing when the clouds overlapping is very slight, 
and in the overlapping regions the momentum spreads of the slow and fast particles are 
small [lO|, i.e., the distribution function roughly is Fa — [S{pz — Po) + 5{j>z +Po)]; where 
Po = 7o^eC and 70 — (l/2)(7max/7min)^^^ — 5. For L waves with = the dispersion 
equation e^z = provides the maximum growth rate T^ax — ^p/("\/27o^^) which is achieved 
at kz^opt = VQ^p/C^VqIo^^), where Vb — c is the particle velocity corresponding to the 
momentum pq. Obliquely propagating L waves {k^ ^ 0) are also generated by the two- 



stream instability. For such waves with kr^ ^ kz the growth rate is ~ ^max III]- The 



characteristic frequency of generated L waves is ~ uopl^jo)^^'^ ~ io^'y)!'^ 

When the amplitudes of generated L waves become large enough linear approximation 
is not valid any more, and the lowest order nonlinear process called quasilinear relaxation 
starts up 



12|. This process may be described by Eq. ([5]), and its characteristic time is 

TQ^^ki^J'^luj,, (14) 

where A is the Coulomb logarithm. The quasilinear relaxation results in that the distribution 
function of particles has a shape of a plateau up to the Lorentz-factor of ~ 70, i.e., is 
~ 7(7"*^ at 7 ^ 7o and nearly zero at 7 > 70. The further generation of A waves is cut off. 

After the stage of quasilinear relaxation the induced scattering of L waves by plasma 
particles becomes the dominant nonlinear process |l2 1 . The characteristic time of this process 
is 

, ^-^^Jf^^ (15) 

where Wl is the energy density of L waves. The induced scattering transfers the wave 
energy from the frequency region where L waves are generated {uj ~ ujp{'~f)^^'^ ^ ^plo^'^ and 
kz ~ kz^opt) to the low frequency region {u ~ a;p(7~^)^/^ ~ ^plo ^^"^ and kz ~ ujp\J (7)/c ~ 
uOp'^y^ /c), i.e., the mean frequency of L waves decreases ~ 70 ~ 5 times. Here, we used that 
(^") ~ 7g for n > and (7"") ~ %^ for n > 1. 
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IV. NONLINEAR DECAY L ^ L' + A 



The process of nonlinear conversion of L waves into A waves, L L' + A (and any other 
three- waves process) is described by the second order current. From Eqs. (jl]) and ([H]) this 
current can be written as 

e^u d 



m dp 

where = + Ef in general case when both L and A waves participate in the process. 

The distribution functions of electrons and positrons created via electromagnetic pro- 
cesses are identical. In such pair plasma all three- waves processes are absent. This is 
because the second order current j''^^ is proportional to e^, and the summation in (fT6ll over 
electrons and positrons gives = 0. However, in the process of outflow of pair plasma 
along the curved magnetic field lines the distribution functions of electrons and positrons 
are shifted with respect to each other to maintain the electric neutrality of the plasma [3] • 
This results in that the mean Lorentz-factors of electrons (7e) and positrons (7p) slightly 
differ, A7 = |7e — 7p| ~ 1. 

For L waves shifted to low frequency region by induced scattering there is no solution 
of the resonant conditions, ujl = ujli + ^"^^ = kii + kL2- Therefore, the process 
L — > Li + L2 is kinematically forbidden, and we can omit all the terms proportional to 
E^{k)E^{k — k') in Eq. f|T6l) . Taking into account that uJi ^ k^c and uja ~ k^c for 
low-frequency L waves and A waves, respectively, from Eqs. (fT3l) and (fT6!) we obtain 



l[UJ 



where 



ujl)E^, = J dkVL'AE^ik - k')E^{k'), (17) 
^ Wn|A7^ 

^L'A V^^J 

is the matrix element of interaction that together with the resonant conditions, 

i^L = ^L' + ^A and = ki' + k^ , (19) 

totally determine the conversion process L L' + A. 

To describe the frequency change in the process L — L' + A it is convenient to introduce 
the ration of the frequency of generated A waves to the mean frequency of the low- frequency 




FIG. 1: Dimensionless frequency of A waves generated in the process L — > L' + ^ as a function of 
z-component of dimensionless wave number of L waves before their decay for different parameters: 
= 0.5, K^' = (sohd hue), = 0.7, K^' = (dashed hne) and = 0.7, K^' = 0.4 
(dash-dotted hne). 

L waves, ^Ia = t^A/("\/2(7^^)^^^t^p) — ooA/iV^lo^^'^i^p)- The results of numerical solution of 
the resonant conditions f|T9l) are shown in Figure 1 where the dimensionless wave number is 



used, K^y = k^^yc/ {y2{'y~^)ujp). We can see that for the typical parameters the frequency 
of generated A waves is about 2 orders of magnitude less then the mean frequency of L 
waves. To clarify our numerical results presented in Figure 1 we can find the approximate 
analytical solution of the resonant conditions (fT9|) for obliquely propagating L waves, Qa ~ 
{K^Y{K^Y /2 ^ 1 at K^' = 0. From this solution it follows that for rather small wave 
numbers of L waves the frequency of generated A waves is always much less than the 
characteristic frequency of L waves. 



Using the standard technique of weak turbulence theory (e.g., l4]), the characteristic 



time scale of nonlinear decay L L' + A for the waves with random phases may be written 



as 



^NLA ~ T/2 , , • (20) 

V. NUMERICAL ESTIMATES AND DISCUSSION 

If we assume that in the plasma frame the frequency of radio emission coincides with the 
mean frequency of L waves generated by the two-stream instability at the distance r^, the 
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characteristic frequency of the radio emission in the pulsar frame is 

VL ~ 10 (5i2M7p7o/Po.i)'/' {Rlnf GHz, (21) 

where Byi is the magnetic field at the neutron star surface in units of 10^^ G, Pq.i is the period 
of the pulsar rotation in units of 0.1 s, i? ^ 10^ cm is the neutron star radius, and M is a 
so-called multiplicity factor that equals to the ratio of the pair plasma density to the density 
of the primary particles j?!. In early papers where only the curvature radiation of primary 
particles was considered as the mechanism of generation of 7-rays the pair multiplicity is 



rather high, M ~ 10^ - 10"^ 



Later, it was shown that generation of 7-rays by inverse 



16j. In this case the value of M 



Compton scattering is essential near at least some pulsars 
may be as small as ~ 1 — 10. For plausible parameters (e.g., 7p ~ 10^, 70 ~ 5, M ~ 10 — 10^, 
5i2 ~ 1, Po.i ~ 1, and TijR = 50) Eq. yields z/^ ~ 2 - 20 GHz that is about an order 
of magnitude larger then the typical frequency of the pulsar radio emission. 

In this paper we discussed the following sequence of processes in the magnetospheres of 
pulsars: (a) development of two-stream instability in collisions of plasma clouds and gener- 
ation of L waves; (b) quasilinear relaxation of two-stream instability; (c) induced scattering 
of generated L waves by plasma particles and decrease of the mean frequency of L waves; 
(d) nonlinear conversion of L waves into low frequency A waves. For plausible parameters of 
the pulsar plasma defined above and for K^^ = 0.3 and Wl ~ 0.1no7o'Tic^ the characteristic 
times of these processes in the plasma frame are = (Fmax)^^ ~ 10^^ s, tqi ~ 10^'' s, 
r/5 ~ 10~^ s, and tjvla ~ 10^^ s, respectively. These times are at least an order of magni- 
tude smaller than the characteristic time of the plasma outflow Tout — ^i/ {cip) — 2 x 10~^ s. 
Therefore, all these processes have enough time to be developed in the magnetospheres of 
pulsars. We have shown that the frequency of A waves is a few ten times smaller than the 
frequency of generated L waves and may be consistent with the observational data on the 
radio emission of pulsars. 

Gedalin, Gruman, and Melrose [l3] proposed another mechanism of pulsar radio emission 
that also involves the two-stream instability and solves successfully the high frequency prob- 
lem at least for the main part of known pulsars. In this mechanism obliquely propagating 
electromagnetic waves are directly generated by the instability. At present, the mechanism 
by Gedalin et al. has an advantage over our mechanism because electromagnetic waves 
generated by their mechanism can escape from the magnetospheres of pulsars (provided 
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cyclotron absorption is unimportant) while it was argued (e.g., [18|) that A waves are sub- 
ject to strong Landau damping and cannot be observed far from the pulsar (but see 19|). 
However, in the simple model used in 18|] it was assumed that the magnetosphere plasma 
is stationary, and the plasma density declines in proportion to without other variations. 
These assumptions are invalid in our case where the pair plasma outflow is nonstationary 



At the edges of 



and strongly non-uniform across and along the magnetic field lines 
the plasma clouds the density falls down sharply and significantly, and A waves may be 
converted into non-damping waves that escape as radio emission. Besides, in the spaces 
between the plasma clouds the density and Landau damping of A waves may be small. A 
study of these effects is beyond the framework of this paper and will be addressed elsewhere. 

In turn, our mechanism has an advantage over the mechanism by Gedalin et al. for pulsars 
with the maximum of the radio spectra at very low frequencies (< 10^ MHz) because the 
mean frequency of radio emission in our mechanism may be an order less than the same in 
the mechanism by Gedalin et. al. 
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